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Introduction to Tidal Streams 
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Abstract Dwarf galaxies that come too close to larger galaxies suffer tidal disrup- 
tion; the differential gravitational force between one side of the galaxy and the other 
serves to rip the stars from the dwarf galaxy so that they instead orbit the larger 
galaxy. This process produces “tidal streams" of stars, which can be found in the 
stellar halo of the Milky Way, as well as in halos of other galaxies. This chapter 
provides a general introduction to tidal streams, including the mechanism through 
which the streams are created, the history of how they were discovered, and the 
observational techniques by which they can be detected. In addition, their use in 
unraveling galaxy formation history and the distribution of dark matter in galaxies 
is discussed, as is the interaction between these dwarf galaxy satellites and the disk 
of the larger galaxy. 


1 An Overview Tidal Streams and Halo Substructure 


The Milky Way is a large, spiral galaxy. It is composed of stars, gas, and dark mat- 
ter within its “halo,” which is the region within which the matter is graviationally 
bound. Like other similar galaxies in the Universe, the Milky Way is also orbited by 
gravitationally bound satellites that contain their own stars, gas, and/or dark matter. 
Globular clusters contain stars, but are generally free of gas and dark matter. There 
are 157 known globular clusters in the Milky Way (Hamid 010) Dwarf galaxies 
contain stars, sometimes gas, and are generally thought to contain large amounts of 
dark matter (though the dark matter content of dwarf galaxies is debated). There are 
about 30 known dwarf galaxies in the Milky Way fawommenm ye Y We gener- 
ally think of globular clusters as one generation of star formation, and dwarf galaxies 
as more massive satellites with ongoing star formation, but the distinction between 
these two classes of objects has become more blurred as we find fainter dwarf galax- 
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ies and more complex globular clusters. Dwarf galaxies can contain globular clusters 
or be orbited by smaller dwarf galaxies. 

When a globular cluster or a dwarf galaxy comes too close to the Milky Way, our 
galaxy's gravitational potential pulls harder on the part of the stellar association that 
is closer to the Galactic center than on the part that is farther away. This differential 
force is called a tidal force, since it is similar to the mechanism that causes the 
ocean tides on Earth. On Earth, it is the gravity of the Moon (and also the Sun) that 
pulls harder on the oceans that face it than on the oceans on the opposite side. This 
difference force stretches the oceans and raises the ocean levels on the sides of the 
Earth that are towards and away from the Moon. 

Tidal forces on stars (or dark matter) in dwarf galaxies and globular clusters pull 
them off off of their satellite of origin, and into tidal streams that can encircle the 
Galaxy. Stars that have been tidally stripped from a satellite of the Milky Way, and 
that are moving more slowly or are closer to the Galactic center than the satellite's 
center of mass, will be pulled into lower energy orbits, which are slightly closer to 
the Galactic center and have shorter orbital periods. Tidally stripped stars that are 
moving more quickly or are farther from the Galactic center than the center of mass 
of the dwarf galaxy or globular cluster will be pulled into higher energy orbits and 
thus orbit more slowly. The physics is analogous to that of planets orbiting the Sun; 
inner planets move faster than outer planets. At the point that the stars are no longer 
trapped by the gravity of the dwarf galaxy (or globular cluster), they orbit the Milky 
Way instead of the center of mass of the satellite. At that point, they are part of the 
tidal tail. As the orbital periods of Milky Way halo satellites are of order hundreds 
of thousands or a billion years, the timescale over which tidal streams evolve is also 
billions of years. 

The content and density of dwarf galaxies is known to vary as a function of 
radius. A typical model for a dwarf galaxy includes ten times more dark matter than 
baryonic matter. The dark matter has a much larger scale radius than the baryonic 
matter, which is more centrally located. There is typically radial structure in the 
stellar populations as well; younger, more metal-rich stars are more centrally located 
than the older stellar populations. These radial gradients in the content of dwarf 
galaxies are imprinted on the structure of tidal streams. 

Stars and dark matter whose populations are distributed at larger radius from 
a dwarf galaxy center will be preferentially stripped from the dwarf galaxy at an 
earlier time (though there is some mixing of material from all radii in the stripping 
process). Material in the tidal streams that is stripped first will generally be located 
farthest from the dwarf galaxy (though it is not a strict time sequence along the 
stream because the stripped stars have a distribution of energies at the time they 
became unbound). Therefore, we expect the dark matter to be stripped first, followed 
by the older populations, followed by the younger populations. This can lead to a 
trend in stellar populations (and chemical abundance) along the stream. 

In addition to radial gradients, dwarf galaxies could have their own possibly com- 
plex or asymmetric shapes; they could be rotating, and often contain globular clus- 
ters within their gravitational potential. Each of these more complex attributes could 
contribute to the structure of the stream. At least some of the globular clusters in the 
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halo of the Milky Way are believed to have been accreted when dwarf galaxies fell in 
and then were assimilated (along with their globular clusters) into the stellar halo. 
Globular clusters themselves, regardless of their origin, may themselves become 
tidally disrupted, though until they become unbound to the dwarf galaxy they would 
be shielded from the tidal forces of the Milky Way. 

The length of the tidal tails depends on the mass and scale radius of the satellite 
(which determine how easy it will be to strip the stars), the path along which the 
satellite orbits the Milky Way, and the length of time the satellite has been orbit- 
ing. While early papers on tidal streams supposed that all tidal stream progenitors 
had been orbiting the Milky Way for the age of the Universe, more recent papers 
recognize that the Milky Way galaxy has likely built up over the age of the Uni- 
verse through merger processes, and the smaller galaxies that are merging now are 
likely to have been captured at later times. Since the tidal forces on a satellite are 
stronger closer to the Galactic center, more material is stripped from the satellite 
and acquired by the tidal tails when the satellite is at perigalacticon. The variable 
rate of mass loss from the progenitor satellite will also lead to structure in the tidal 
tails. Models for generating the tidal stream from the Sagittarius dwarf galaxy, for 
example, often segregate the tidal stream stars by the perigalactic passage during 
which they were stripped. 

Observed tidal streams of stars could have a known progenitor which still sur- 
vives, despite the fact that some of its stars have been removed. Streams might also 
have no known progenitor either because the satellite from which it formed has been 
completely disrupted into the tidal stream, or because it has yet to be discovered (for 
example it might be hidden behind the Galactic disk). Stars that have been stripped 
by tidal forces are subsequently incorporated into the stellar halo, while still main- 
taining some kinematical and chemical memory of their origin. 

Before the discovery of tidal streams, the stellar halo was modeled as a smooth 
distribution of stars, fit by a power law profile (p ~ r-™, where p is the density, r is 
the distance from the Galactic center, and à is a number close to 3.5; the distribution 
was often additionally “squashed,” in the sense that the density fell of more rapidly 
in the direction perpendicular to the Galactic plane). This was sensible when it was 
expected that halo stars were formed during a rapid collapse of a monolithic gas 
could from which the Milky Way was formed 1962). It even 
seemed sensible when the subsequent model of (SZ 1978), which 
described the buildup of stars in the halo through mergers, was proposed; it was 
assumed that there were so many individual tidal streams that formed the halo that 
the overall density was relatively smooth. 

We know now that tidal streams cause significant density substructure in galactic 
halos. Most of the earlier measurements of the parameters in the original power law 
models depended on which lumps of tidal debris were probed by the observations. 
While it is now clear that at least some of the halo built up through recent mergers, a 
portion of the halo stars could be part of a smooth density distribution. If there was 
an early period of rapid merging, the result would still be a smoothly varying spa- 
tial distribution of stars that also appears well-mixed in velocity. It is not currently 
known whether all of the stars in the halo were born in other, smaller galaxies that 
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later merged with the Milky Way, or whether some of the stars were formed during 
an initial, rapid collapse of the Milky Way, as proposed by ELS in 1962. Power law 
models are still fit to the shape of the stellar halo, but a conscious choice must be 
made about whether to fit the density as a smooth component plus significant over- 
densities in large tidal streams, or whether to average the halo density over large 
volumes and thus include the large tidal streams in the fit. 

In fact, the majority of the tidal debris that has so far been discovered is between 
20 and 50 kpc of the Galactic center. It may be that interior to the Solar Circle 
the stellar halo consists of so many tidal streams that have been dispersed over so 
many orbits that the result is indistinguishable from a rapid collapse. Farther than 50 
kpc, the satellites are subject to very weak tidal forces, very long orbital times, and 
may be very late to fall into the Milky Way's gravitational well. Therefore, these 
very distant satellites might not have formed tidal tails. Studies indicate that the 
fraction of halo stars in substructures is a few tens of percent within 30 kpc of the 
Galactic center, and the fraction increases at larger radii crt Hu and 
references therein) 

The Milky Way stellar halo may have a smooth, well-mixed component. But 
it also contains dwarf galaxies, globular clusters, tidal streams, and "clouds." The 
"clouds" are density enhancements in the halo that are many kiloparsecs across and 
could not possibly be gravitationally bound, but are also not drawn out into linear 
structures like tidal streams. These structures are thought to be the result of a minor 
mergers that come in on highly eccentric orbits, which take them near (or through) 
the center of the Milky Way. Simulations show that these mergers could result in a 
pile-up of stars at the apogalacticon of the orbit, which would look like a large cloud 
when viewed from inside the Milky Way. 

In addition to the Milky Way, tidal streams have also been observed in other 
galaxies (see, e.g., Chapters 8&9 of this volume). These streams are most notable 
in large galaxies close to the Milky Way, such as the Andromeda galaxy (M31; dis- 
cussed in Chapter 8). Note that we are using the term “tidal stream" to mean the stars 
pulled from a satellite during a minor merger with a much larger galaxy. This is a 
distinct phenomenon from the ejection of stars or gas that may occur during a major 
merger between two galaxies of similar size. Major mergers can form structures that 
give a stream-like appearance. A classic example of this latter type of structure is 
the the Antennae Galaxies (NGC 4038 and NGC 4039). 


2 The Discovery of Tidal Streams in the Milky Way 
2.1 Dwarf galaxies and moving groups 


During the last two decades of the 20" century, following the landmark paper of 
Searle & Zinn| (1978), the Galactic structure community began to consider the pos- 


sibility that dwarf galaxy mergers could be shaping the morphology of the Galac- 
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tic halo. [Rodgers et al] Observed metal-rich, main sequence A stars near the 
south Galactic pole. Because these stars were far from the star-forming regions of 
the Galactic disk, and main sequence A stars have a limited liftime (and therefore 
cannot be found far from their place of birth), they concluded that the A stars could 
have been the result of an encounter between the Milky Way and a dwarf galaxy. AI- 
though twenty years later these stars were identified as more likely blue stragglers in 
the thick disk (Girard et al. (2004), the paper was important in pushing astronomers 
to look for halo moving groups, that would result from dwarf galaxy mergers. At 
that time, blue straggler stars were thought to exist only in globular clusters, where 
the high stellar density would result in stellar mergers and/or larger binary fraction. 

It was quite reasonable in the 1980's to assume that all main sequence A stars 
must be young. In fact, a decade later, Preston et al. (1994) were still operating under 
the assumption that blue stragglers were preferentially found in the high stellar den- 
sity environments of globular clusters; field main sequence stars were still strongly 
argued to be young. We now believe that some (or all) blue stragglers are the result 
of mergers of binary stars. This process produces a high mass main sequence star 
long after the initial formation of the stellar population. This process can occur at 
any time, in any stellar population, regardless of stellar density. 

Though evidence for moving groups of halo stars continued to grow through the 
late 1980's and early 1990’s, the first solid evidence for tidal streams in the Galactic 
halo was presented by|Majewski et al] (1996). This paper is the third in a series that 
measures the proper motions and radial velocities of stars near the north Galactic 
pole, within 8 kpc of the Galactic plane, and finds a retrograde moving group of 
20 stars with a net velocity down towards the Galactic plane. In addition to this 
dominant group, they find evidence for additional velocity substructure. Though this 
particular moving group has not yet been connected to other tidal substructure, many 
observations have since shown the same phenomenon; if a small volume of halo 
stars are studied, significant structure in the velocities of the stars is observed. At the 
same time, theoretical studies of the growth of galaxies in cold dark matter (CDM) 
cosmology concluded that galaxies grew through mergers 1988). 

The discovery of the Sagittarius dwarf galaxy, which appeared to be in the pro- 
cess of tidal disruption [1994), added to the evidence for hierarchical 
merging. Here was a modern-day example of a dwarf galaxy in the process of being 
assimilated into our galaxy. Milky Way globular clusters were identified as possi- 
bly stripped from the Sagittarius dwarf spheroidal AEAEE ood. 
Ard|Lvaden-Bell & Lynden-Bell suggested that other halo globular clusters 
could represent “ghostly streams” from dwarf galaxy mergers with the Milky Way, 
though they were not able to to show that any of their suggested streams were “un- 
equivocally real.” 

In the last decade of the 20" century, people were talking about kinematic discov- 
ery of tidal streams in the halo, but there was no expectation that the tidal streams 
would be detectable in stellar density. There were expected to be a large number 
of tidal streams, that came in from random directions, and each contributed a very 
small fraction of the density at any position in the halo. The tidal debris from dwarf 
galaxies would be detectable as a moving group in the Galactic halo for more than 
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a gigayear (Johnston et al} |1995) and would remain on great circles for more than 
a gigayear (Johnston et al. 1996). Because models of dwarf galaxy tidal disruption 
showed that the Sagittarius dwarf galaxy could not have survived on its present or- 
bit for the age of the Milky Way, models with mass-to-light ratios of 100 were put 
forward to explain it's survival until the present day 1997). Later, it 
would be suggested that dwarf galaxies could be deflected into their present orbits 
by other dwarf galaxies or from the large number of dark subhalos expected to be 
present in the Milky Way’s halo, so it was not required for the Sagittarius dwarf in 
particular to have survived on its current orbit for a substantial fraction of the age of 
the Universe. 

It should be noted here there there is a vast literature on galaxy mergers in the 
second half of the twentieth century, and it is a bit difficult to define when that 
community started to recognize halo substructures that are similar to the Sagittarius 
dwarf tidal stream in our galaxy. The term “tidal stream” was used to describe the 
distortions induced in galaxies during major mergers, rather than the way we are 
using it in this book: to describe a group of stars that have been torn off of a much 
smaller satellite during an accretion event. Conference proceedings showing very 
faint tidal features in external galaxies in the late 1990’s, but we particularly high- 
ligh a paper by [Shang et al] (1998), that was published at about the same time the 
Sagittarius dwarf and its tidal stream were being discovered in the Milky Way. In 
this paper they discovered a dwarf galaxy in NGC 5907 this is of similar size to the 
Sagittarius dwarf galaxy. They suggested that this dwarf galaxy could be causing 
the previously observed warp in H I. They also identified a faint, elliptical ring that 
appeared to be the remains of a dwarf galaxy that had been tidally disrupted. 

successfully identified kinematic substructure in the stellar 
halo near the Sun, where accurate 3D velocities were available. They examined 
the 3D kinematics of 97 low metallicity giant and RR Lyrae stars within 1 kpc of 
the Sun, and found that about 1096 of them had similar velocities. We show a re- 
reduction of Hipparcos catalog data used in that paper in Figure [I] Helmi et al. 
concluded that 1096 of all metal-poor halo stars outside of the solar circle were 
accreted from one dwarf galaxy merger with the Milky Way, under the assumption 
that the local halo represented the halo as a whole. While the moving group of stars 
was confirmed by (2009), the idea that it represented a component 
of the entire halo was not. Subsequent papers refer to the moving group as a tidal 
stream from a satellite that was accreted 6-9 billion years ago, and which comprises 


about 5% of the local halo (Kepley et al.,/2007). 


This success in kinematic detection of halo substructure helped to launch the 
"Spaghetti Survey" which aimed to find kinematic substruc- 
tures through the Milky Way halo using high latitude pencil beam surveys; small 
fields were imaged in colors that allowed the photometric identification of red giant 
and blue horizontal branch stars, for which spectra were later obtained. The final 
sample included 101 well-observed halo giant stars, of 
which 20% were known to be members of tidal streams. The survey concluded that 
between 10% and 100% of halo stars are in (presumably accreted) substructure. 
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Fig. 1 These panels show the local velocity distribution and angular momentum (Jz and Jperp = 


Jg --Jg). in Galactic XYZ coordinates, of nearby stars in the Hipparcos survey, following 


(1999). Here, X points in the direction from the Sun to the Galactic center, Y is in 
the direction of solar motion, and the coordinate system is right-handed. The sign of the angular 
momentum opposite to the original paper, which used a left-handed coordinate system. The blue 
filled circles identify stars in a coherent halo moving group. 


2.2 Halo substructure in the era of large surveys 


The Sloan Digital Sky Survey (SDSS; York et al. 2000) was designed to map the 
large scale structure of the Universe by measuring the positions of galaxies. The dis- 
tances to about one million of the brightest galaxies in 10,000 square degrees would 
be measured from spectroscopically determined redshifts. The brightest galaxies 
would be selected from a high accuracy imaging survey of of the sky in five optical 
passbands: u, g,7,i, and z. Although the extragalactic program provided the primary 
science driver for the survey, the project provided tremendous opportunities for all 
other areas of astronomy, and was particularly influential in fostering the field of 
tidal streams. 

First light for the SDSS was obtained in May of 1998, when the imaging survey 
began with observations of a 2.5?-wide stripe of data along the Celestial Equator 
(since observations on the equator did not require the telescope to track). Unex- 
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Fig. 2 This plot shows how low surface gravity BHB stars can be separated from high surface 


ravity BS stars using ugr colors. The colors of model stars, with varying metallicity, are from 
(1998). Also shown are colors of BHB and BS stars in Pal 5, and several other stars 


with spectroscopically determined surface gravities. The heavy black line is the empirical curve 
that roughly separates BHBs and BSs by color. Figure 10 off Yay et all (2000). 


pected density substructure in the Galactic halo was apparent from the very begin- 
ning, though because these were early data reductions it took us a some time to 
believe that what we were seeing was not an observational artifact. 

The first (gravitationally unbound) halo substructure in density was found from 
blue horizontal branch and blue straggler stars in the first 2.5?-wide stripe of 
SDSS imaging data on the Celestial Equator.{Yanny et al] selected blue stars 
(—0.3 < g — r < 0.0), that were initially expected to consist entireley of blue hor- 
izontal branch stars; main sequence A stars and blue stragglers were generally 
thought to be absent in field halo populations. This paper showed that there were 
two large over-densities of stars - one in the north Galactic cap, and one in the south. 
It was also shown that the overdensities included both blue horizontal branch stars 
and blue stragglers. The blue stragglers were two magnitudes fainter, and could be 
roughly separated from the BHB population with ugr colors (Figures[2]and[3). The 
structure in the north Galactic cap was also discovered in SDSS RR Lyrae stars at 
about the same time 2000). Because these substructures were diffuse 
and unlikely to gravitationally bound, we assumed they were two separate, tidally 
disrupting structures in the Galactic halo. 

It turned out, though, that these two apparently separate structures are two cross 
sections through the Sagittarius dwarf tidal stream, which is the most prominent 
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Fig. 3 These three wedge plots show the blue straggler (BS) and blue horizontal branch (BHB) 
stars along the Celestial Equator, as selected from early data from the SDSS. Galactic coordi- 
nates are in square brackets and Equatorial coordinates are in parentheses; the solid line shows 
the intersection of the Celestial Equator and the Galactic plane. The top plot shows all stars with 
dereddened colors 0.8 < u* — g* < 1.5 and —0.3 < g* — r* < 0.0; the lower plots show BHBs (left) 
and BSs (right), selected by color as shown in Figure[2] The denser arcs at 195? < œ < 230° and 
20? < a < 40° show the positions where the Sagittarius dwarf tidal stream crosses through the 
Celestial Equator. The fainter arcs (larger g magnitude) are the intrisically fainter BS stars, and the 
same distance as the BHB stars, which are about two magnitudes brighter. Notice that the BHB 
stars in the southern (upper left) part of the Sagittarius stream appear to be offset in angle from the 
BS stars. It turns out that the Cetus Polar Stream (CPS [2009), which is very rich in 
BHB stars, happens to cross through the Celestial equator at almost the same position as the Sgr 


dwarf tidal stream, and most of the southern BHBs that we originally thought were part of the Sgr 
stream are actually from the CPS. Figures 3, 11, and 12 tom anm stall bog. 


tidal stream in the Milky Way. This was shown by (2001b), who traced 


the tidal stream across the sky with carbon stars. Two years later, this tidal stream 
was traced across the sky in M giant stars identified photometrically in the Two 
Micron All-Sky Survey (Figure[4] Majewski et al. 2003). The discovery of a large 
tidal stream in the Milky Way galaxy was quickly followed up by the discovery of 
a large tidal stream in the Andromeda Mn [baia etal] (2001a), the Milky Way's 
nearest neighbor. 
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Fig. 4 The positions of M giant stars within 10° of the Sgr dwarf tidal stream, selected from 
2MASS, are shown. The (Xs,,. Gc. Ysg;Gc) coordinates are the Sagittarius dwarf stream coordinates, 
where the stream is in the plane Zsggc = 0. Because the Sgr stream is close to polar, this plot is 
roughly equivalent to Galactic Z vs. X, where the center of the Galaxy is at (0,0), the Sun is at 
X — —8.5, Y is the direction of solar motion, and Z points towards the north Galactic pole. The 
Sgr dwarf is shown by a large filled circle and a line indicating the direction of its motion. The M 


giants trace the path of the Sgr dwarf tidal stream all the way around the Galaxy. Figure 11 from 
2003). 


Thus was born the field of tidal streams. Although the Sagittarius dwarf galaxy 
tidal stream remains the most prominent and well-studied tidal stream in the Milky 
Way, even fifteen years later we are still struggling to model the tidal disruption of 
this dwarf galaxy. At first, many researchers assumed that the tidal stream from the 
Sagittarius dwarf galaxy was an anomoly - that this was the only tidal stream that 
would be identified as a stellar overdensity rather than a co-moving group of stars. 
The discovery of density substructure was, after all, unexpected. 

It was not long, however, before additional substructure began to be discovered. 
found an overdensity of RR Lyrae stars in the Virgo constellation. 


Tidal tails were found extending tens of degrees from the Pal 5 globular cluster 
(Odenkirchen et al], 2001 2002). Previous to these papers, globular 


cluster “tidal tails” referred to stars outside the tidal radius of the globular cluster, 
that were identified as “wings” in the density profiles (Grillmair et al} 1995]. 
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In 2002, Newberg et al. made the seminal discovery that the density of main 
sequence turnoff stars in the Galactic halo is quite lumpy, and is a poor fit to a 
smooth, power law density profile. The density of turnoff stars on the Celestial 
Equator (Figure[2.2) includes a wealth of halo substructure. In this slice through the 
Celestial Equator, one can see the Sagittarius dwarf tidal stream in both the north- 
ern and southern Galactic hemispheres; the Pal 5 globular cluster, which turns out 
to have long tidal tails; the overdensity in Virgo that was identified by [Vivas et al] 
(2001): stellar concentrations towards the Galactic center that could be due to a 
smooth component of the stellar spheroid, but also include stars in what was later 
named the Hercules-Aquila Cloud (2007); a controversial over- 
density near the anticenter commonly referred to as the “Monoceros Ring;” a more 
distant overdensity south of the Galactic plane called the “Triangulum-Andromeda 
Stream” (identified in the original paper as Monoceros in the south); and two con- 
centrations of closer stars near the anticenter (one in the north and one in the south) 
that appeared be due to the disk but at the time did not fit any disk model. 

Also in 2002, Ferguson et al. showed that there was spatial and metallicity sub- 
structure in the halo and outer disk of M31 (Figure[6), using red giant branch stars 
in a panoramic photometric survey of this galaxy. M31 is the closest large, spiral 
galaxy to the Milky Way, and has provided a complementary view of tidal structure 
in spiral galaxies. While the Milky Way allows us to see stream structure in three 
dimensions, it is much easier to see the overall halo structure in external galaxies. 

The discovery of lumpy outer halos containing tidally disrtupting dwarf galaxies 
in both the Milky Way and Andromeda (M31) bolstered the idea that galaxy forma- 
tion happens through a series of mergers. The lumps in the halo are evidence that 
small mergers continue to happen at late times, and the Sagittarius dwarf galaxy is 
an example of a merger that is still happening today. Since these discoveries, many 
more tidal streams have been discovered in both galaxies, and also in galaxies out- 
side of the local group. 

An amazing map of the density of turnoff stars, color-coded by distance (see 
Figure[7), was published by (2009); it shows the most prominent 
tidal streams in the north Galactic cap, spread out across the sky. In this image one 
can see not only the Sagittarius dwarf tidal stream, but also a fainter tidal stream 
next to it that may also be associated with Sagittarius, but so far no model has been 
able to simultaneously fit both of the “bifurcated” streams. The Virgo Overdensity 
is also conspicuous, as is the “Monoceros Ring,” and another smaller tidal stream 
called the “Orphan Stream,” because it has no known dwarf galaxy progenitor. The 
Orphan Stream progenitor either has yet to be discovered or has been completely 
disrupted. 

Also apparent, nearer to the Galactic center, is the Hercules-Aquila “cloud.” 
Overdensities that are not long and thin like tidal streams, but instead look like 
large, amorphous clumps are referred to as “clouds.” These could be produced by 
tidal debris from a dwarf galaxy on an orbit that takes it near the center of the Milky 
Way galaxy (Hemauist & Ouind (198. The tidal disruption is severe as the dwarf 
galaxy goes through the center of the Milky Way, where the tidal forces are strong. 
The stars that have been tidally stripped then orbit back and forth through the Galac- 
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tic center, spending a larger fraction of their time at apogalacticon, where they are 
observed as amorphous "clouds" of stars. Note, however, that it has been suggested 
that at least one of these clouds (the Virgo Overdensity) could be the disrupted re- 
mains of a dwarf galaxy that is on a highly eccentric orbit, but just recently passed 
perigalacticon and thus does not consist of debris piled up at the apocenter of its 
orbit (Carlin et all 2012). 

While this section gives an historical account of the major observational mile- 
stones in the beginning of the field of tidal streams, the succeeding chapters of this 
book describe in detail the discoveries that have since been made. When spatial sub- 
structure in the halo was first discovered, Don York, from the University of Chicago, 
said we had discovered “Eggen’s spaghetti sky.” The same Olin Eggen that is fa- 


mous for a paper suggesting the Milky Way formed from the rapid gravitational 
collapse of a cloud of gas is also famous for seeing clumps of 
stars in imaging data, and in particular introducing the concept of “moving groups” 
of stars [loss [Ee oen & Sandaad [1959) - stars that are in a similar vol- 
ume of space and moving in the same direction. Although his moving groups are 
likely from tidal disruption of open clusters in the disk rather than dwarf galax- 
ies and globular clusters the halo, the image of a sky filled with stars from tidally 
disrupted clusters mixed together like a bowl of spaghetti is a common analogy. 
make a similar analogy to a “can of worms.” However, the 
width of these stellar streams is astonishing. The tidal stream from the Sagittarius 
dwarf galaxy, for example, has a width of about 6 kpc. For reference, the distance 
from the Sun to the Galactic center is about 8 kpc, and the radius of the disk is about 
15 kpc. The width of 6 kpc is a substantial fraction of the width of the Milky Way. 
Although tidal streams from globular clusters might create “spaghetti” in the sky, 
dwarf galaxies like Sagittarius are made of much thicker pasta! Perhaps we should 


call it the “rigatoni sky!” This is of course written in jest, but one should not forget 
the enormous spatial extent of some of these streams and clouds. 


3 Observational Techniques for Finding Substructure 


Tidal streams can either be identified as an overdensity of stars in the halo, or by 
finding stars with similar locations and velocities (using angular momentum allevi- 
ates the need for co-location). The latter technique will allow the identification of 
tidal streams with lower surface density, a longer time since disruption of the satel- 
lite, or which are closer to the Galactic center where the stellar density is higher. 


3.1 Standard Candles and Photometric Parallax 


The most straightforward way to identify tidal streams by stellar density is to de- 
termine the three dimensional positions of halo stars and then look for localized re- 
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gions that have densities higher than their surroundings. Measuring positions in 3D 
requires that one know the distance to each star, so particular stars with known ab- 
solute magnitude (standard candles) must be selected. Often, these standard candles 
are selected photometrically based on their colors, but they might also be selected 
from spectroscopy or as variable stars. 

Figure [8] shows a color-magnitude diagram (CMD) of stars from the Palomar 5 
globular cluster, showing typical features for an old, metal-poor stellar population. 
The brighter stars (which allow us to probe the distant halo) include blue horizontal 
branch stars (BHBs), RR Lyrae variables, red giant branch stars, M giants, blue 
stragglers, K giant stars, and turnoff stars (note that not all of these types are present 
in Pal 5!). There are advantages and disadvantages to using each type of star. 

Horizontal branch stars are good standard candles because they occupy a very 
narrow range of absolute magnitude. Blue horizontal branch stars (BHBs) can be 
preferentially selected by blue color, though with some contamination from blue 
straggler stars which have the surface gravities of main sequence stars. However, 
BHBs are only present in very old, metal-poor stellar populations, and there are 
relatively few stars of this type. Red horizontal branch stars are more difficult to 
identify due to confusion with giant branch stars and much fainter main sequence 
stars. 

RR Lyrae stars are even better standard candles, but are harder to identify (usually 
through multi-epoch variability studies), and also appear only in intermediate-age, 
fairly metal-poor populations. RR Lyrae stars are even sparser tracers of the stellar 
density than horizontal branch stars. Additional observations (spectra and/or time- 
series photometry) is required to determine the distance to each RR Lyrae star. 

M giant stars are extremely bright, but only found in metal-rich stellar pop- 
ulations, which are rare in the halo of the Milky Way. In fact the only verified 
tidal stream that has been traced with M giant stars is the largest and most ap- 
parent overdensity: the Sagittarius dwarf tidal stream (see Chapter 2). The Mono- 
ceros Ring, Triangulum-Andromeda cloud (see Chapter 3), and Pisces overdensity 
have also been studied in M giants. 

K giant stars are a natural choice for tracing stellar densities because they are 
more numerous than other stars of similar intrinsic brightness, and present in all 
old stellar populations. The main drawback to using K giant stars is that spectra 
(or specialized filters that are particularly sensitive to the gravity of K stars; see, 
oe. Majewskict all are required to separate K giants from the much more 
numerous K dwarf stars in the disk, which have similar apparent magnitudes as K 
giant stars in the halo. Additionally, the absolute magnitude depends on the age and 
metallicity, and is a strong function of color (red giant branches are nearly vertical 
in the CMD); the measured distances will depend heavily on the model used to 
calibrate them. 

Blue stragglers are best separated from blue horizontal branch stars using spectra, 
but they can be roughly separated if there is data in the u, g, and r filters. These 
stars are not usually used to trace stellar populations because they are not good 
standard candles (absolute magnitude is a strong function of color), they are fairly 
rare, and they are not as bright as giant stars so they do not probe as far in distance. 
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However, they are often observed in tidal streams, two magnitudes fainter than the 
blue horizontal branch stars, and must be considered particularly when using BHB 
stars to trace stellar density. 

Turnoff stars are much more numerous than any of the other types of stars con- 
sidered above, but they are by far the worst standard candles - spanning at least two 
magnitudes at about constant color. When using this type of star to measure density 
along the line-of-sight, it is advantageous to use statistical photometric parallax, 
which is described in Section 3.3. 


3.2 Matched Filter Techniques 


Matched filter techniques are used in many applications to optimally extract signal 
from noise. In general, one convolves the observed, noisy signal (as a function of 
either space or time) with a template (filter) that looks like the expected signal. This 
will enhance signal and surpress noise. 

In this case, we would like to find the angular density of stars in a tidal stream 
population (that have a known distribution in color and apparent magnitude), while 
filtering out stars that are from other disk or halo populations (for which it is also 
possible to construct an expected distribution in color and apparent magnitude). 

Matched filtering was first used to detect tidal streams by Reckoai call (2002), 
who discovered the tidal tails of the Pal 5 globular cluster in SDSS data. The 
matched filter technique described here follows the description in that paper. In 
this version of the matched filter, the cluster stars that are in regions of the color- 
magnitude diagram where few background stars are expected are weighted more 
heavily; a minimum variance estimate (defined below) for the number of cluster 
stars is derived, rather than a maximum signal-to-noise estimate. By more heavily 
weighting the stream detections in regions of low background, systematics from 
poor knowledge of the background is reduced. 

One defines a normalized probability in color-magnitude space that describes 
where one expects to find stars in the stellar population of the stream (including 
primarily turnoff, giant, and blue horizontal branch) as N;;,. In the case of Pal 5, the 
CMD of the GC itself, normalized so the total probability integrated over the whole 
color-magnitude diagram is one, is used for Ny. In addition, one defines a normal- 
ized probability in color-magnitude space for the background, Npg. In general, this 
background depends on Galactic latitude and longitude. Then, an estimate for the 
number of stream stars, œ, in a given solid angle can be found by minimizing the 
following expression for 77: 


, 


" IN(color, mag) — |&Ns; (color, mag) + Nos (color, mag)] }* 
pg (color, mag) 


color,mag 


where N(color,mag) is the number of stars in a particular color, magnitude bin. 
The sums are over all of the bins that have a non-zero number of stars expected 
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in the stream. The background must be constructed so that the expected number of 
stars is also non-zero in all of these bins. This equation can be solved for the œ 
that corresponds to the minimum Y? by setting the first derivative with respect to œ 
equal to zero: 


Y N(color,mag)Ns;(color,mag) 
Npg(color,mag) 


— Na, (color, mag) 


color,mag 


N2. (color,mag) 
color,mag Npg(eolor,mag) 
In the limit of infinitely small bins, the first term in the numerator can be replaced 
by a sum of N,;, (color, mag) / Np, (color, mag) for each of the observed stars in a 
particular solid angle. The second term in the numerator and the denominator can 
be replaced by an integral over color-magnitude space rather than a sum over the 
individual bins. 

In the case of Pal 5 tidal tails, Ns; was derived from stars with high membership 
probability in the globular cluster. The background distribution, N54, was derived 
from stars in several square degree patches of actual SDSS data near the globular 
cluster. Because the background distribution was derived from the actual survey 
data, it includes features due to both the survey itself (including increased color 
errors and decreased completeness near the limiting magnitude of the survey), and 
the density distribution of stellar populations in the Milky Way. The tidal tails of 
Pal 5 were expected to be at about the same distance from the Sun as the cluster 
itself across the region surveyed, so no shift in apparent magnitudes of stream stars 
was considered. Because the parts of the background color-magnitude diagram that 
varied across the part of the sky surveyed were at colors and magnitudes where there 
were few stream stars, a varying background distribution was not considered in this 
paper. Rockosi et al. constructed a smoothed density plot of the tidal tails of Pal 5 
with a resolution of 0.09? by finding the optimal value of œ as a function of sky 
position. 

Grillmair and Johnson (2006) modified this matched filter technique so that it 
could be applied to the discovery of tidal streams when the distance to the streams 
and the stellar population is not known a priori. The description here follows Grill- 
mair (2009). 

The probability distribution for the stream stars is generated from globular clus- 
ter stars observed in SDSS data. The mean g — r and r — i colors are determined for 
a fixed set of g magnitudes in each of the clusters; colors at other magnitudes are 
interpolated. The distribution can be shifted brighter or fainter to account for differ- 
ent distances to the tidal debris. The relative number of stars expected in the stream 
as a function of apparent magnitude (the luminosity function) is determined from 
example globular clusters with very deep photometry available, and modified by the 
expected completeness as a function of apparent magnitude for the SDSS survey. It 
is broadened in color by the amount of the expected photometric error as a function 
of apparent magnitude. An additional broadening with o — 0.02 magnitudes is also 
applied to account for the intrinsic dispersion in the colors of giant branch stars. 
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This new technique makes it possible to create a set of 2D pictures of the sky, 
that show the matched-filter density maps of different stellar populations, at different 
distances. These can then be searched for linear structures that indicate the presence 
of a tidal stream. One creates “filters” that consist of the expected color-magnitude 
distribution of stream stars in the SDSS, at different distances and using different 
globular cluster distributions as the template. A variety of background distributions, 
distances, and stellar populations, are used to create the 2D pictures. One can scan 
through a distance series of these pictures, for a fixed stellar population, to pick out 
linear structures that might vary in distance as a function of position in the sky. 


3.3 Statistical Photometric Parallax 


If one observes objects whose stellar type can be identified from photometry, and 
which therefore have a known absolute magnitude, then distances to these objects 
can be determined by photometric parallax (see section 3.1). An example where 
photometric parallax for a large number of stars was used to measure the density 
of stars in the Milky Way is [furiéetal] (2008). They estimated the distances to 
48 million SDSS stars of type G and later, under the assumption that they were 
main sequence stars (which most of them undoubtedly are), and using a photometric 
calibration from (2008). From the estimates of the individual distances, 
the positions of each star are derived, and from the resulting stellar densities the 
structure of the Milky Way (including tidal substructures) was derived. 

In contrast, statistical photometric parallax (Newberg, uses a statistical 
approach to determine the underlying spatial density of stars from the distribution 
of observed apparent magnitudes. It differs from photometric parallax because dis- 
tances are not determined individually for each star. This technique has been used 
to determine the spatial distribution of stars in the Sagittarius dwarf tidal stream, 
using turnoff stars from the SDSS (Newby et a1][2013). Although turnoff stars have 
a wide variety of absolute magnitudes, and are thus not very good standard candles, 
they are quite useful for density mapping because they are much more numerous 
than giant stars, and can be seen at larger distances than fainter, lower mass main 
sequence stars. Statistical photometric parallax is most successful for measuring 
structures with comparable (or higher) stellar density than the background at their 
location. 

The basic concept of statistical photometric parallax is demonstrated in Fig- 
ure [2.2] There is a narrow, radial feature at RA=229° at an apparent magnitude of 
g ^» 21. This feature is the globular cluster Pal 5, whose stars are all essentially the 
same distance from the Sun, but the apparent magnitudes of the turnoff stars in the 
figure span an apparent magnitude range from g ~ 20 to g ~ 22.5. The large range of 
apparent magnitude is a result of the large range of absolute magnitudes of turnoff 
stars. In fact, one could make a histogram of the apparent magnitudes of the Pal 5 
turnoff stars, shift it by the distance modulus to Pal 5, and obtain the distribution of 
absolute magnitudes of color-selected Pal 5 turnoff stars. 
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(2011) showed that the distribution of absolute magnitudes of color- 


selected turnoff stars were similar for all halo globular clusters observed in the 
SDSS. This surprising observation is apparently a result of the age-metallicity rela- 
tion for the Milky Way galaxy. Older stellar populations should have fainter, redder 
turnoffs because a larger number of the more massive stars have had time to evolve 
away from the main sequence. On the other hand, the older stellar populations are 
typically more metal-poor, and metal-poor main sequence stars are brighter and 
bluer than metal-rich stars of the same mass. It is surprising, but apparently true, 
that in the halo of the Milky Way these two effects nearly cancel so that the dis- 
tribution of absolute magnitudes of blue turnoff stars is essentially the same for all 
halo stellar populations. Although the color of the turnoff will be slightly different 
for different ages and metallicities, the absolute magnitude distributions are surpris- 
ingly similar. 

This age-metallicity conspiracy provides a great opportunity to apply statisti- 
cal photometric parallax to measure the spatial density of halo stars, including the 
substructure from tidal streams. The observed halo substructure, for example as 
depicted in Figure [2.2] represents the actual density structure, convolved with the 
absolute magnitude distribution of halo turnoff stars, and then modified by the ob- 
servational constraints. 

To determine the spatial density, we first construct a parameterized spatial den- 
sity model. Given a particular set of parameters, one constructs the expected density 
distribution over a particular volume, derives from this the expected apparent mag- 
nitude distribution as a function of angular position in the sky (taking into account 
the known absolute magnitude distribution of turnoff stars), and then applies obser- 
vational biases such as completeness. This expected distribution of stars can then be 
compared with the observed distribution of stars. The parameters can be optimized 
to most closely match the observed data. For example, used a 
maximum likelihood technique that was optimized with conjugate gradient descent 
to measure the density of tidal streams, and in particular the Sagittarius dwarf tidal 
stream, in the halo. 

As with matched filtering, proper application of this technique includes many 
subtleties. For example, if turnoff stars are selected in a narrow color range, then the 
distribution of stars selected will depend on the color errors. Typically color errors 
increase considerably near the observation limits, so that stars with a wider range 
of colors are sampled at the faint end of the survey. Because there are few stars 
bluer than the turnoff, and many intrinsically fainter main sequence stars redder 
than the turnoff, the distribution of absolute magnitudes of apparently faint stars in 
a magnitude-limited survey includes a larger number of intrinsically faint stars. 

A benefit of the statistical photometric parallax technique is that once the spatial 
density of a halo substructure is determined, it is possible to statistically remove 
stars with the spatial distribution of one halo structure or substructure from the sam- 
ple, so that the remainder of the sample can be studied. Note, however, that it is not 
possible to determine which stars are actually a part of that structure or substructure. 
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4 Co-Moving Groups of Stars 


Several halo substructures have been initially identified as group of stars of simi- 
lar type (RR Lyrae, BHB, etc.), in a small volume, that have similar velocities. In 
some cases, these co-moving stars have been associated with a stellar overdensity 
that is clearly the result of tidal stripping. In other cases the spatial extent of the re- 
lated stars is unknown. In this case the association may be referred to as an Element 
of Cold Halo Substructure (ECHOS, [2009). The disadvantage to 
this technique is that it requires that the stars have known radial velocities as mea- 
sured from stellar spectra, or measured proper motions, or both. The advantage is 
that as an increasing amount of astrometric and spectroscopic data become available, 
this technique will allow us to detect and characterize halo substructure with much 
lower density contrast that result from small mergers, old mergers, mergers where 
the satellite has been widely spatially dispersed (for example by passing through the 
center of the Milky Way), or tidal streams in high density regions of the Milky Way. 
Kinematically identified halo moving groups are covered in more detail in Chapter 
2: 

Halo tidal streams have not been identified inside the Solar Circle (within 8 kpc 
of the center of the Milky Way). This could be because these more central halo stars 
formed in situ or because they are due to a large number of well-mixed mergers. If 
the central regions formed from mergers, it might be possible to separate the stars 
into groups with similar angular momentum. The angular momentum coherence 
will persist for much longer than spatial coherence, and all of the stars that fall in 
on a particular satellite have similar angular momentum (see Chapter 6). 


5 Chemical Tagging 


A star’s chemical composition is determined by the chemical composition of the 
gas cloud from which it formed. Stars in the same open cluster or globular cluster 
are thought to be coeval and therefore have similar chemical composition (though 
at least some globular clusters contain more than one stellar population). If one as- 
sumes that most stars are born in star clusters, and that stars maintain a memory of 
the gas cloud in which they were formed because their outer layers contain the same 
chemical ratios as those gas clouds, then one could attempt to group stars by the 
cluster in which they were born (Freeman & Bland-Hawthornl boos. This requires 
that one make very accurate measurements of a few tens of elemental abundances 
with high resolution spectroscopy. The hope is that the chemical abundances of pri- 
mordial clusters are different enough that they are distinguishable in the observable 
abundance space. Chemical tagging, particularly of the Galactic disk, is the primary 
science driver of the GALAH survey (De Silva et al. 2015). 

In contrast to star clusters, dwarf galaxies generally contain more than one 
episode of star formation. While studies of the angular momentum of halo stars 
will elucidate the history of satellite mergers that built up the halo, studies of stel- 
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lar chemistry can tell us about the individual star clusters from which the stellar 
populations are composed. In particular, the lowest metallicity stars tell us about 
the conditions in the early Universe, when these stars formed, and about the earliest 
chemical enrichment processes. Intermediate metallicity stars require more than one 
episode of star formation (which is not guaranteed if gas is ejected from low mass 
dwarf galaxies in the first epoch of star formation), and teach us about the environ- 
ments in which the stars were formed throughout the evolution of the Universe. It 
would be incredible if, at some time in the future, we could build up a record of both 
the accretion event in which each halo star became a part of the Milky Way, and the 
star formation event in which that star was born, within the accreted satellite. 

Currently, only the simplest chemical tagging techniques have been applied to 
the detection of tidal streams (for example, the detection of debris from the globular 
cluster Omega Cen via its unique abundance signatures by Pefaiowski et all Boi, 
though streams are often characterized by their chemical composition. 

Of particular interest is the variation of chemical composition as a function of the 
distance from the progenitor. Because the outer (higher velocity dispersion) layers 
of a satellite are tidally stripped first, the tidal streams can tell us about the density 
structure of the progenitor dwarf galaxy. We expect that higher metallicity stars are 
are more centrally located in the dwarf galaxy, and therefore will be stripped at later 
times and be preferentially located nearer to the progenitor dwarf galaxy in the tidal 
stream. In fact, metallicity gradients have been reported in the Sagittarius dwarf tidal 
stream (see Chapter 2). 


6 Constraining Dark Matter, the Formation of the Milky Way, 
and Cosmology with Tidal Streams 


6.1 Rapid collapse vs. Hierarchical Mergers 


In the late 20" century, there was rapid evolution in our understanding of galaxy for- 
mation. Eggen, Lynden-bell & Sandage are credited with putting forward the mono- 
lithic “rapid collapse” picture of the formation of the Milky Way, in which halo stars 
were formed during the collapse of a proto-Galactic gas cloud. This picture was well 
summarized in a review by [Eogen et all (1973): “The oldest Population II objects 
(subdwarfs, globular clusters, metal-poor RR Lyrae stars) formed when the galaxy 
was extended and far from equilibrium, before much metal enrichment occurred. 
The remainder of the galaxy continued contraction to a disc-like configuration, ap- 
proximately in centrifugal equilibrium, and dissipated the necessary energy before 
forming stars; rapid metal enrichment took place. This is the disc population. The 
residual gas settled to a very flat disc in which the star formation is still going on; 
this is the Population I." 


Searle & Zinn (1978) painted a competing picture of galaxy formation, in which 
they proposed "that the gas from which the clusters and stars of the outer halo 
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formed continued to fall into the Galaxy for some time after the collapse of its cen- 
tral regions had been completed; that the interactions of the infalling gas dissipated 
much of its kinetic energy and gave rise to transient high-density regions in which 
the halo stars and clusters formed; that these regions dispersed even while they 
underwent chemical evolution; that the stars and clusters that had formed within 
them eventually fell into dynamical equilibrium with the Galaxy and constitute its 
present outer halo, while the gas lost from these protogalactic star-forming regions 
was eventually swept into the galactic disk." In other words, the halo stars are stars 
formed in other, smaller galaxies and clusters, which later merged with a galaxy that 
had already formed. This grew into the "hierarchical merger" scenario, in which the 
Milky Way grew to its current size through mergers of smaller galaxies and stellar 
associations. 

The hierarchical merger picture has been borne out both by simulations (e.g. 
and by the fact that we are observing 
recent and present-day accretion of small galaxies into the Milky Way halo. Note, 
however, that the majority of the merging may have happened at early times, and 
thus the rapid collapse model is not completely incorrect. Also, hydrodynamical 
simulations suggest that galactic disks may have formed through steady flow of 
cold gas into massive dark matter halos (Dekel et al. [2009). 

Because many stars currently in the Milky Way could have been formed in 
smaller galaxies and later accreted, the star formation history in our galaxy may 
be different from the history of structure formation. For example, the composition 
of stars that were formed outside the Milky Way does not tell us anything about the 
composition or density of the gas in the Milky Way at the time that star was formed. 
There could be an epoch of star formation that is followed by an epoch of merging. 


6.2 Constraining Dark Matter with Tidal Streams 


In the standard ACMD model (Spergel et al] |2015) of the Universe, most of the 


mass in galaxies like the Milky Way is made of of dark matter, that interacts with 
itself and the rest of the mass in the Universe only through gravity. Around Milky 
Way-sized dark matter halos, there could be thousands of dark matter subhalos that 
are the size of typical dwarf galaxies. Because this matter does not interact with 
light, we cannot see it directly. Instead we infer its existence from the motions of 
objects (such as stars and galaxies) that we can see, and from gravitational lensing 
of objects such as galaxies and QSOs. 

Although the stars in the Milky Way halo (often called spheroid stars) repre- 
sent only 1% of the stars in the halo, they provide an archaeological record of our 
galaxy's merger history. A significant fraction (or maybe all) spheroid stars were 
formed in other, smaller, galaxies that fell into and were ripped apart by the Milky 
Way. Any dark matter associated with those galaxies is incorporated into the dark 
matter halo of the Milky Way. At least some of the Milky Way globular clusters 


Introduction to Tidal Streams 21 


are believed to have come into the Milky Way with dwarf galaxies that are now 
disrupted. 

The disruption of dwarf galaxies into particular Milky Way tidal streams is typ- 
ically modeled using N-body simulations. In these simulations, the Milky Way is 
often assumed to be a fixed potential, and the progenitor satellite is modeled as a 
symmetric distribution of N (say 10,000 or 100,000) particles of equal mass, whose 
positions and velocities then evolve through time due to gravitational forces from 
each other and from the fixed Milky Way potential. Since satellites at 20 to 50 kpc 
from the Galactic center take of order a billion years to orbit once, these simulations 
are typically run for billions of years. However, they cannot be trusted further back 
in time than about 3 billion years, since the gravitational field of the Milky Way 
galaxy should evolve with time - particularly at early times in the Universe. These 
N-body simulations can be used to constrain the properties of the progenitor, the 
orbit of the satellite from which it was derived, and the gravitational potential of the 
Milky Way itself. 

Alternatively, tidal streams and halo substructure can be modeled in the cosmo- 
logical context using N-body or hydrodynamic simulations of galaxy formation in 
general. In this case, the galaxy(ies) and tidal streams will be produced, but will not 
match the observed galaxies and tidal streams we observe in the Universe today, but 
they can be used to ask whether the types of structures we see are consistent with 
the models. Examples of large simulations of galaxy formation include: the Aquar- 


ius Project (Springel et al. (2008), Via Lactea (Diemand et al, |2008) and Illustris 
(Vogelsberger et al. 


Because the stars in the tidal streams were ripped from their progenitor satel- 
lites by the Milky Way's gravity, and their subsequent orbits are determined by the 
Milky Way's gravity, tidal streams hold great promise for teaching us about the grav- 
itational potential of the Milky Way. Because the majority of the mass in the Milky 
Way is believed to be make of dark matter particles, tidal streams could tell us the 
density distribution of dark matter. 

Although in principle one could determine the three dimensional position and 
velocity of any star in the Milky Way, stars in tidal streams are the only ones that 
we know where they were in the past. We know that tidal stream stars were all once 
co-located on the same progenitor satellite. This is the property of tidal stream stars 
that makes them potentially very valuable tracers of the Galactic potential. 

Dark matter in the Milky Way galaxy can be described by its overall shape. In 
addition to this overall (smoothly varying) density distribution, the radial velocities 
of the stars in dwarf galaxies suggest that these satellites contain a large fraction 
of dark matter eS ios d Moreover, cosmological models of structure 
formation produce galaxies with many more “subhalos” than there are observed 
dwarf galaxy satellites; there could be many dwarf galaxy sized lumps of dark mat- 
ter, that do not contain stars, in orbit around the Milky Way. 

Tidal streams have been used to constrain dark matter in all of these contexts (see 
Chapter 7). The Sagittarius tidal stream in particular (see Chapter 2) has been used 
to infer a triaxial shape for the Milky Way's dark matter halo. Dark matter subhalos 
can weakly interact with stream stars creating a wider stellar stream, or strongly 
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interact with stream stars, throwing them out of the stream and creating gaps. The 
width (or alternatively velocity dispersion) of tidal streams gives some indication of 
the mass of the progenitor, which can in principle be compared with the number of 
stars observed to determine the amount of dark matter in a progenitor dwarf galaxy. 


7 Disk Response to Tidal Interactions 


When large dwarf galaxies (or dark matter subhalos) pass close to the Milky Way, 
the disk can also respond to the gravity of the dwarf galaxy 
[2009). Historically, this mechanism was studied as the ori- 
gin of the observed warp in the Galacitic disk, and in particular the warp was thought 
to be the result of an interaction with the Large Magellanic Cloud (LMC). After the 
discovery of the Monoceros Ring (See Chapter 3), disk response to passing satellites 
was proposed as a mechanism to explain the presence of this highly controversial 
structure in the plane of the Milky Way 
Poog). Later work focussed on the Sagittarius dwarf galaxy, which is known to be 
close to passing through the disk, as the most likely progenitor (Purcell et al) [201 1). 

Surveys of stars have corroborated this picture of disk oscillations. Coherent 
disk velocity substructure was observed in the SDSS (Widrow et al. [2012), the RA- 
dial Velocity Experiment (RAVE;|Williams et al], 2013) , and the Large Area Multi- 
Object fiber Spectroscopic Telescope (LAMOST; borg. Oscillations 
in the disk density have also been observed as a function of height above the plane 
[Yanny & Gardned DOL. and distance from the Galactic center bord. 
The observed density waves are in reasonable agreement with the simulations of 
(2013), who calculated the affects of the Sagittarius dwarf galaxy on 
the disk. 

Evidence for wavelike disk oscillations has also been found in the Andromeda 
galaxy. Andromeda’s NE clump (see Chapter 8), 40 kpc from the galaxy center, 
is thought to be composed of disk stars, possibly tossed away from the disk plane 


by an encounter with the dwarf galaxy progenitor of Andromeda’s Giant Stellar 
Stream 2015). It is unclear how the corru- 
gations in Galactic gas and star forming regions observed in our galax (Quiroga, 
ona. and in others Matthews & Usonl 
2008), are related to the disk oscillations that are being observed in Milky Way 
disk stars, but the subhalos have also been proposed to explain the ripples in the 


gas (Chakrabarti & Blitz 2009). One wonders whether, if disk oscillations could be 


caused by galactic subhalos, if they also induce spiral structure. 
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8 Future Prospects 


A very dynamic picture of present-day galaxies is emerging. The stellar halos of 
large galaxies appear to be composed (in part or in total) of stars that have been 
tidally stripped from smaller galaxies that have fallen into the larger galaxy’s gravi- 
tational field. Massive subhalos (including dark matter-only subhalos, if they exist) 
could be exciting waves in galactic disks. We expect that subhalos could produce 
a time-varying gravitational field in the Milky Way, that could among other things 
produce gaps in tidal streams. There is even a possibility that very massive infalling 
dwarf galaxies could be significantly moving the center of mass of the larger galaxy, 
affecting not only the disk but also the paths of the dwarf galaxies and associated 
tidal streams (Gómez et al, [2013]. 

As ongoing and future surveys provide us with large numbers of radial veloc- 
ities and proper motions of Galactic stars (the largest of these surveys will come 
from the ESO Gaia satellite), it will be possible to map a larger fraction of halo 
substructure kinematically (see Chapter 5), and to better constrain the properties of 
the streams we have already identified. This better data will allow us to trace tidal 
stream substructure with higher fidelity, and thus help clarify which stars in tidal 
streams originated on the same Milky Way satellite. We will be able to map the 
spatial extent of kinematically identified streams (which have density contrast too 
low to follow by density alone), and the kinematic substructure of tidal streams that 
have been identified by density contrast. 

The Sagittarius dwarf tidal stream (see Chapter 2) is presents illuminating ex- 
ample of how far we have come in understanding tidal streams, and how far we 
still have to go. This stream is the largest and best studied of the Milky Way tidal 
streams, but continues to defy a clear explanation of how its structure evolved to 
match these observations. We have traced the tidal stream more than 360? around 
the Milky Way in K/M giant, BHB, and turnoff stars. We have measured the veloc- 
ities of stream stars and their velocity dispersion, mapped the stream width and the 
density profile along the stream, measured the chemical compositions and chemical 
gradients along the stream, compared the properties of the stars in the stream with 
those in the progenitor galaxy, and identified apparent bifurcations of the stream into 
two pieces in both the north and south Galactic caps. However, we are still wonder- 
ing what combination of Galactic potential and dwarf galaxy progenitor would give 
rise to these observations. On the observational side, it is not a certainty that all of 
the identified parts of the stream are associated with a single dwarf galaxy progeni- 
tor. On the theoretical side, we struggle with how spatially and temporally complex 
the model for the Galactic potential should be. We also note that the properties of the 
progenitor dwarf, including rotation, substructure, satellites, and the radial profiles 
of the stellar and dwark matter constituents, can give rise to observable complexity 
in the resulting tidal stream. On the one hand, tidal streams have the potential to 
provide us with a wealth of information about our galaxy and the satellites that have 
fallen into it. On the other hand, they require us to understand a large number of 
complexities to reach definite conclusions. 
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Better understanding of tidal streams will lead to better understanding of galaxy 
evolution, and to better constraints on the distribution of dark matter, both in the 
progenitor dwarf galaxies that are thought to have very high mass-to-light ratios, 
and in the Milky Way galaxy itself. Techniques for constraining the spatial struc- 
ture of dark matter are currently in development (See Chapters 6 & 7). We still face 
challenges in understanding the dynamical systems that we are modeling, and in 
parameterizing the dark matter distribution, that could be quite complex both spa- 
tially and in time. However, the wealth of streams also provide a large number of 
constraints that one day might unravel the dark matter puzzle. 
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Fig. 5 Wedge plot (g vs. RA) showing the density of 334,066 SDSS turnoff stars (0.1 < (g—r) < 
0.3) on the Celestial equator. This illustrates the discovery that the Milky Way's stellar halo 
contains a tremendous amount of density substructure. The extent and nature of this substruc- 
ture remains an active area of research, but much of the observed substructure arises from dwarf 
galaxies that have fallen into the Milky Way and been tidally disrupted. The angle shows right 
ascension, and apparent g* magnitude increases along radius. The solid line indicates the inter- 
section between the Celestial Equator and the Galactic plane. The Pal 5 globular cluster is evi- 
dent at (a, g) = (2295,21); the radial extent in the diagram is due to the spread in the absolute 
magnitudes of turnoff stars. The Sagittarius stream is at (/, b, g) = (341°,57°,22.5) in the north 
Galactic cap, and at (/,b,g) = (167°, —54°,21.5) in the south Galactic cap. The Virgo Overden- 
sity is at (/,b,g) = (297°, 63°, 20.0). The Hercules-Aquila Cloud juts out from the overdensity 
at (/, b, g) = (52°, —32°, 20.4). There is a region of high extinction near œ = 60° that affects the 
identification of turnoff stars; the radial feature there is not an actual density substructure. Near 
the Galactic anticenter, 70° < œ < 135^, overdensities can be seen north of the plane at g = 15 
(attributed to the disk in the original paper), south of he plane at g — 17 (identified as the Mono- 
ceros ring in the south bylibats er all 2003), north of the plane at g = 19.4 (the Monoceros Ring), 
and south of the plane at g = 20 (the Triangulum-Andromeda stream, identified as the Monoceros 
ring in the south in the original paper). The feature labeled (/,b, g) = (6°,41°,20) is the only 
overdensity that could plausibly represent the smooth distribution of spheroid stars that had been 
previously expected. 
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Fig. 6 We show the surface density of redder red giant branch (RGB) stars in Andromeda (M31). 
The outer ellipse has a semimajor axis length of 55 kpc; the inner ellipse has a semimajor axis 
of 27 kpc and is well outside of the region containing the optical disk of M31. Halo substructure, 
including the giant stellar stream at (&, n) = (0.5, — 1.5) and lumpy structure surrounding the disk. 
This image showed that the halo of Andromeda contained * nificant density substructure, just as 


is seen in the Milky Way. Figure 2 from[Ferguson et al] 12002). 
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Fig. 7 This spectacular image of density substructure in the Milky Way's stellar halo is known as 
the "Field of Streams." Declination vs. right ascension (increasing towards the left) is shown for 
SDSS turnoff stars in the north Galactic cap. This is an RGB composite where blue represents the 
brighter (closer than about 15 kpc) stars, and red represents the faintest (more distant than about 
25 kpc) stars. Note that there appear to be two branches to the Sagittarius Stream, and a multitude 
of halo substructures. The circles indicate the positions of new Milky Way satellites discovered in 
SDSS data; two are faint globular clusters and the others are faint dwarf galaxies. Image credit: V. 
Belokurov and the Sloan Digital Sky Survey. 
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Fig. 8 The color-magnitude 
diagram shows stars that can 
be used to trace Galactic 
structure. The black points 
are SDSS stars within four 
arcminutes of the Pal 5 glob- 
ular cluster, and represent the 
distribution of stars in the halo 
population. In a magnitude- 
limited survey, brighter stars 
sample the galaxy to larger 
distances, but they are not 
very numerous. Turnoff stars 
have been popular for trac- 
ing halo structure because 
they are easily identified by 
color, and are the brightest of 
the more numerous main se- 
quence stars. The drawback to 
turnoff stars is that they have a 
large spread in absolute mag- 
nitude, so individually they 
are not good standard candles. 
This difficulty is somewhat 
mitigated by the fact that, for 
halo stars in the Milky Way, 
the absolute magnitude distri- 
bution is not a strong function 
of population. Therefore, sta- 
tistical photometric parallax 
can be used to recover the un- 
derlying density distribution 
from the positions and ap- 
parent magnitudes of turnoff 


stars. Figure 1 of 
(2013). 
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